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Abstract 

When solitary waves are characterized as homoclinic orbits of a finite-dimensional 
Hamiltonian system, they have an integer- valued topological invariant, the Maslov 
index. We are interested in developing a robust numerical algorithm to compute 
the Maslov index, to understand its properties, and to study the implications for 
the stability of solitary waves. The algorithms reported here are developed in the 
exterior algebra representation, which leads to a robust and fast algorithm with 
some novel properties. We use two different representations for the Maslov index, 
one based on an intersection index and one based on approximating the homoclinic 
orbit by a sequence of periodic orbits. New results on the Maslov index for solitary 
wave solutions of reaction-diffusion equations, the fifth-order Korteweg-De Vries 
equation, and the longwave-shortwave resonance equations are presented. Part 1 
considers the case of four-dimensional phase space, and Part 2 considers the case of 
2n— dimensional phase space with n > 2. 



1 Introduction 

Hamiltonian evolution equations in one space dimension, such as the nonlinear Schrodinger 
(NLS) equation, fifth-order Korteweg-De Vries (KdV) equation, longwave- shortwave reso- 
nance (LW-SW) equations, have the property that their steady part is a finite-dimensional 
Hamiltonian system. For such systems, solitary wave solutions can be characterized as 
homoclinic orbits of the Hamiltonian ordinary differential equation (ODE). The spectral 
problem associated with the linearization about a given homoclinic orbit, in the time- 
dependent equations, then leads to a parameter-dependent family of linear Hamiltonian 
systems. The advantage of these Hamiltonian structures is that the linear and nonlinear 
Hamiltonian systems have global geometric properties that aid in proving existence of the 
basic solitary wave and in understanding its stability as a solution of the time-dependent 
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equation. Our interest in this paper is in a particular geometric invariant - the Maslov 
index of homoclinic orbits. 

The study of the stability of solitary waves using the Maslov index was pioneered in 
the papers by Jones [27J and Bose & Jones [5J. The linear stability of standing wave 
solutions of a spatially-dependent NLS equation is studied in [57]. The linearization about 
a steady solution results in a linear A— dependent Hamiltonian system with two degrees 
of freedom of the form (\1.2\\ and A a spectral parameter. Geometric methods are then 
used to determine the Maslov index, and it is used to prove an instability result. Gradient 
parabolic PDEs of the form 



u t = d\u xx + f u (u, v) 
v t = d 2 v xx + f v {u,v) 



'1.1) 



are considered in [6J, where d\ and d 2 are positive parameters, f(u,v) is a given smooth 
function with gradient (f u ,fv)- Linearizing about a steady solution (u(x),v(x)), and in- 
troducing a spectral parameter leads to a coupled pair of linear second-order ODEs which 
can be put into the standard form (11.21) . with the asymptotic property (11. 5ft and A the 
spectral parameter. Since the PDE is a gradient system it is sufficient to restrict the 
spectral parameter to be real. Singular perturbation methods are then used to determine 
the Maslov index, which in turn is related to stability. A key feature of this work is the 
analysis of the induced system on the exterior algebra space A 2 (^ 4 )- 

Many of the most interesting solitary waves are only known numerically and therefore a 
numerical approach to the Maslov index is of interest. It is the aim of this paper to develop 
a numerical framework for computing the Maslov index of homoclinic orbits. Once the 
solitary wave solution is known, analytically or numerically, it is the linearization about 
that solitary wave which encodes the Maslov index. Therefore, the starting point for 
developing the theory is the following class of parameter-dependent Hamiltonian systems 



Ju x = B(x, A)u, uel 4 , xeM, Ag 
where J is the standard symplectic operator on IR 4 

-I 

1 

and B(x, A) is a symmetric matrix depending smoothly on x and A. Let 

A(x,X) = J _1 B(a;,A) . 



1.2) 



(1.3) 



;i-4) 



The fact that A(x, A) is obtained from the linearization about a solitary wave suggests 
the following asymptotic property. It is assumed throughout the paper that 



Aoo(A) = lim A(x, A) 



;i.5) 
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and that Aoo(A) is strictly hyperbolic for an open set of A values that includes 0. 

The theory applies to linear Hamiltonian systems on a phase space of dimension 2n with 
n any finite natural number. In Part 1 attention is restricted to the case of 4— dimensional 
phase space which simplifies formulae, and general aspects of the case n > 2 are given in 
Part 2 [IB]. 

The Maslov index is a winding number associated with paths of solutions of (11.21) . in 
particular, paths of Lagrangian planes. A Lagrangian plane in IR 4 is a 2— dimensional 
subspace of M 4 , say span{zi,z 2 }, satisfying (Jzi,z 2 ) = 0, where (•, •) is a standard inner 
product on M 4 . 

Suppose A is fixed and on the interval a < x < b consider a path of Lagrangian planes 

[a, b] h-> Z(x, A) = [ Zl (x, A) | z 2 (x, A)] G M 4x2 , 

satisfying Z x = A(x, A)Z for a < x < b. The Maslov index of this path is a count of the 
number of times this path of Lagrangian planes has a non-trivial intersection with a fixed 
reference Lagrangian plane. A precise definition is given in $21 

A byproduct of the present theory is some new observations about the properties of the 
Maslov index, which in turn are useful in computation. In the numerics, the exterior 
algebra formulation is also advantageous. We give formulas for different representations 
of the Maslov index for Lagrangian planes on A 2 (^ 4 ) (and for any n in Part 2 [18]), and 
present a general algorithm that works - in principle - for any dimension n. However, the 
dimension of A n (^ 2n ) increases rapidly with n and so the algorithm is most effective for 
low dimensional systems. The algorithm is constructed so that the manifold of Lagrangian 
planes is attracting. Numerical results are presented in this paper for the cases of M 2 and 
M 4 . 

In order to develop a numerical framework for the Maslov index, one of the first difficulties 
is defining the Maslov index. Although it is easy to give a rough definition, making it 
precise depends greatly on the context and a surprising number of special difficulties and 
cases arise. 

We will appeal to two constructions of the Maslov index. The first is based on an in- 
tersection index between the Lagrangian path and a reference plane. This definition was 
used in Maslov's original work, and was developed further by Arnold [3] and Duister- 
maat [25]. It is this definition that is used by Jones [27] and Bose & Jones [6], taking 
the Lagrangian path to be a path of unstable subspaces and taking the reference plane to 
be the stable subspace at infinity. 

Independently, Chen & Hu [19] give two constructions of the Maslov index of a ho- 
moclinic orbit. Their first definition is based on an intersection index and is equivalent 
to the definition in [2"Tfo] . Their second definition is based on a Fredholm index of (11.21) 
viewed as an operator in a function space on the real line. However, this latter definition, 
although equivalent to the definition based on an intersection index, is not convenient in 
the numerics since it requires to compute a spectral projector of an infinite dimensional 
operator and then to determine a Fredholm index. 
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In this paper, the definition of the Maslov index based on an intersection index is extended 
by introducing an explicit and computable formula for the crossing form. This theory is 
developed in §H1 

Our second method for computing the Maslov index is to approximate the homoclinic orbit 
by a sequence of periodic orbits, apply the Maslov index for periodic orbits, and then take 
limits. There does not appear to be any loss of generality in using periodic approximates. 
Vanderbauwhede & Fiedler [S] prove that homoclinic orbits in Hamiltonian systems 
(as well as reversible systems) can be approximated as the limit of a sequence of periodic 
orbits. 

The Maslov index for periodic orbits has been widely developed because of its interest 
in semi-classical quantization (e.g. |26ll21f34f42ll39U38] and references therein). In [TB] a 
new numerical scheme is developed to compute the Maslov index of hyperbolic periodic 
orbits, and Chardard [13] proves under suitable hypotheses that if the periodic orbit is 
asymptotic to a homoclinic orbit, the Maslov index converges to the Maslov index of the 
limiting homoclinic orbit. This approach ties in with existing schemes for computing the 
basic solitary wave, where the solitary wave is approximated by a periodic orbit and then 
a spectral method is used for computation. 

The computational framework for the Maslov index is illustrated by application to four 
examples. The first is a tutorial example on IR 2 , where the details can be given explicitly. 
It is a scalar-reaction diffusion equation with an explicit localized solution. The second 
example is a coupled reaction-diffusion equation which also has an explicit solution. The 
third example is solitary wave solutions of KdV5. The fourth example is the LW-SW 
resonance equation which arises in fluid mechanics and consists of a NLS equation coupled 
to a KdV equation. This latter example has two new interesting features: the spectral 
problem is on a six-dimensional phase space, and for appropriate parameter values it has 
a Maslov index which is a non-monotone function of A. 



2 Linear Hamiltonian systems and Lagrangian subspaces 

A Lagrangian subspace can be represented by a Lagrangian frame: a 4 x 2 matrix of rank 
2 




where X and Y are 2x2 matrices satisfying 

Y T X = X r Y . (2.2) 

There is a correspondence between elements of the unitary group U(2) and Lagrangian 
subspaces. When X T X + Y T Y = I then X ± iY are unitary matrices. The determinant 
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of a unitary matrix lies on the unit circle. This property suggests defining the angle 



_ det[X-iY] 
6 ~det[X + iY]- {Z - 6) 

Along a path of Lagrangian subspaces, this angle will change, and the winding of this 
angle is the basis of the Maslov index. 

Let Z(x), a < x < b be any smooth path of Lagrangian subspaces. If the path is a loop: 
Z(6) = Z(a), then there is an integer associated with the path: the number of times the 
induced path on the unit circle, represented by e 1K ^, encircles the origin. Define the angle 
associated with the path by 



i<x) = det\K(x)-iY(x)] 
det[X(x) + iY(z)] ' 

Then the Maslov index of the path is the integer 



(2.4) 



hi(b) — K,(a) 

Maslov(Z, k) := — , (2.5) 

where k here is viewed as the lift from S 1 to R. 

Geometrically Span(Z) — > Maslov(Z, k) induces an isomorphism between the homotopy 
group of A(n) at a point A and the group of integers (Z, +). 

The Maslov index can also be defined for an arbitrary path of Lagrangian subspaces by 
introducing the idea of an intersection form. This approach to the Maslov index goes back 
to Maslov and Arnold [3] . The key to the definition in this case is the use of an inter- 
section form or crossing form. Here we will follow the construction in ROBBINS [4*2lf4l3] and 
Robbin & Salamon jlQ]. Modulo a choice of orientation these definitions are equivalent. 

Let V be a fixed reference plane, represented by a Lagrangian frame. For example a typical 
choice for the reference plane is 

(, 

V = 

In the case of homoclinic orbits a natural choice for the reference plane is the stable or 
unstable manifold at some value of x. 

Suppose, for some value of x, denoted xo, there is a simple intersection between the 
reference plane V and the path: that is, Z(xo) H V is one dimensional. The intersection 
index at xq is determined by the sign of T, the crossing form, defined by 

T(Z, V, x ) = (3Z'(x )P, Z(x )(3) vol . (2.6) 

Here and throughout the paper vol is taken to be the standard volume form on M 4 , 

vol = ei A e 2 A e 3 A e 4 . 
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In the formula (12 .6j) . 

Z'(*„) : 



dx 



Z{x), 

-xo 

and (3 E MP 1 , where MP 1 is the one-dimensional real projective spac^Zl The parameter 
(3 determines the linear combination of the columns of Z(xq) which span the intersection 
subspace 

Z(x ) H V = span{£} , £ := /3 lZl + /3 2 z 2 , 

where Zi,Z2 are the columns of Z{xq). At each simple intersection between and the 
path Z{x) the sign of the intersection form is ±1. Adding the intersection indices over the 
path gives the Maslov index 

Maslov(Z, V) = J2 signr(Z,V,x ). (2.7) 

a<xo<b 

This formula assumes that Z(a)C\V = {0} and Z(b)C\V = {0} and that there is only a finite 
number of Xq where Z(x ) H V ^ {0}. When the endpoints have non-trivial intersection 
the formula can be modified to contribute a half-integer for each end intersection (see 
page 831 of [10]). Non- regular intersections are not generic in the one-parameter family, 
and so can be eliminated by perturbing A. 

Geometrically, two paths image(Z), image(W) in A(n) with the same endpoints and such 
that Maslov(Z, V) and Maslov(W, V) are well-defined, are in homotopy if and only if there 
exists V such that Maslov(Z, V) = Maslov(W, V). 



3 The Evans function associated to the set of linear Hamiltonian systems 

In order to compare the number of eigenvalues of (11.21) with the Maslov index, we will 
use the Evans function to determine eigenvalues based on the setup in Alexander, 
Gardner & Jones [TJ, adapted to the symplectic setting in Bridges & Derks [8], 
restricted to the case of M 4 . 

Consider the linear system of ODEs, 

u z = A(x,A)u, uel 4 , (3.1) 

where A(x, A) = J _1 B(a;, A) and B(x, A) is symmetric and depends smoothly on x and 
A. In general A can be complex but in this paper it will be restricted to be real. Assume 
that A(x, A) tends exponentially fast to a matrix Aoo(A) when x — > ±oo. 



i 

MP 1 = { p e R 2 : (3 + 0, /?~c/3, c G R \ {0} } . 
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Define the stable and unstable subspaces of Aoo(A) by 

EL, (A) := {u G K 4 : lim e Aoo(A)x u = 0} 

x— >+oo 

and 

E^A) :={u6l 4 : Jmi^ e Aoo(A)x u = 0} 

E S (A) (E^j(A)) is the direct sum of the generalized eigenspaces associated with the eigen- 
values of Aoo(A) with negative (positive) real part. The matrix Aoo(A) is said to be hyper- 
bolic if 1R 4 = E^(A) © E^ 3 (A); equivalently if Aoo(A) has no purely imaginary eigenvalues. 
Purely imaginary eigenvalues of A oc (A) are associated with the essential spectrum. The 
essential spectrum is 

c"ess = {A G C | Aoo(A) is not hyperbolic} 

(3.2) 

= { A G C : det[A oc (A) - \kI\ = for some k G R } . 

We will assume throughout that A ^ a css . Then the Hamiltonian symmetry of Aoo(A) 
gives that dimE^(A) = dimE^(A) = 2. 

Let A 2 (R 4 ) be the vector space of 2— vectors in IR 4 . There is an induced system from (13.11) 

\] x = A (2) (x,A)U, UgA 2 (K 4 )- (3.3) 

Let cr + (A) be the sum of the eigenvalues of Aoo(A) with positive real part, and let cr_ (A) be 
the sum of the eigenvalues with negative real part. Then there are solutions U + (x, A) and 
U~(x, A) of (13.31) with maximal decay as x goes to — oo and +oo respectively satisfying 

lim^_ 0O e- CT +W a; U + (x, A) = C + (A) G A 2 (» 4 ) , (3.4) 

and 

\im x ^ +00 e-°-^ x U-(x, A) = C(A) G A 2 (K 4 ) , (3.5) 
where C ± (^) are eigenvectors 

A^ ) (A)C ± (A) = ( r ± (A)C ± (A). (3.6) 

The eigenvalues cr±(A) are analytic functions of A and so the eigenvectors C ± (A) can be 
chosen to be analytic as well. 

A value A G K \ o~ ess is called an eigenvalue if the stable solutions U~(x, A) and unstable 
solutions U + (x, A) have nontrivial intersection. Eigenvalues are detected by the Evans 
function [I] which is defined by 

D{\) vol = U-(z, A) A U + (x, A) G A 4 (K 4 ) • (3.7) 
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The Evans function is independent of x and is an analytic function of A [T]. Analyticity 
assures that the zeros of D(X) are isolated. In the definition 03 .7p . the property of (11.21) . 

Trace(A(x, A)) = , (3.8) 

has been used. This property follows since A = J _1 B with J skew-symmetric and B 
symmetric. 

Because of the Hamiltonian structure, the Evans function is invariant under exponential 
scaling of the following form. Let 

\J ± (x, A) = e~ a± ^ x U ± (x, A) , 

Then the scaled functions satisfy 

\% = [A®(x,\)-*±(\)1\Xj ± , 

but the Evans function becomes 

D(X) vol = e (CT - (A)+a+(A))x, tn(x, A) A U + (x, A) = U~(x, A) A V + (x, A) , 

since cr_(A) + cr+(A) = Trace(A OC) (A)) = 0. 

An explicit expression for the entries of A^ as a function of the entries of A, when 
n = 2, is given in Appendix El It is natural to ask whether there is an induced symplectic 
structure on A 2 (M 4 ). For example, can the induced system be written in the form 

j( 2 )U, = B( 2 )(a;,A)U, 

where J^ 2 ** is the induced matrix from J on A 2 (IR 4 )- However, this is not the case. The 
most significant obstacle is the fact that is not invertible. The precise relation between 
j(2) ) g(2) anc [ _^ (2) j g gj ven j n Appendix El 

4 U ± (a;, A) represent paths of Lagrangian planes 

The paths of stable and unstable subspaces U ± (x, A) (or their scaled versions U ± (x, A)) 
are paths of Lagrangian subspaces. 

Let $(x, s, A) be a fundamental solution matrix for (II. 2p . that is, 

J$ x . = B(x,A)$, $(s,s,A) = I, 
and $(x, s, A)$(s, t, A) = t,X). Define the stable and unstable subspaces for each 

x e R [II], 

E s (x , A) = {u G M 4 : lim $(x, x )u = 0} = span |col(U"(x , A))} , 
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and 

E u , A) = {u G R 4 : lim $(x, x )u = 0} = span (col(U + (a;o, A))} . 
Both subspaces define invariant vector bundles over R. This means that 

E s (x,A) = $(x,s,A)E s (s,A) and E"(x, A) = s, A)E u (s, A) . (4.1) 

Moreover 

lim E s (x, A) = E S (X) and lim E u (x, A) = E" (A) . 

x — >+oo X — >— oo 

If E n (x, A) is Lagrangian for some x then it is Lagrangian for all x. This observation is 
implicit in (14.11) but a direct proof can be given as follows. When u, v G E u (x, A) are 
solutions of (11.21) . 

-^(Ju, v) = (Ju XJ v) + (Ju, v x ) = (B(x, A)u, v) - (u, B(x, A)v) = , 

using symmetry of B(x, A). Hence the value of (Ju, v) is an invariant of (II. 2p for any pair 
of vectors u, v G R 4 : 

(Ju(x, A), v(x, A)) = (Ju(x , A), v(x , A)) , Vx . 
But u, v G E u (x, A) and so lim^^oo u(x, A) = lim ;c _ > _ 00 v(x, A) = 0, therefore: 
lim (Ju(x, A), v(x, A)) = =^ (Ju(xo, A), v(xo, A)) = . 

X — > — oo 

Therefore, E u (x, A) is a Lagrangian subspace for any x. A similar proof confirms the result 
for E s (x, A). Another proof is to use Montaldi's Theorem [37] on Lagrangian planes in R 4 
and a sketch is given in Appendix O 



5 An example in R 2 

Before proceeding to the full definition and properties of the Maslov index for paths of 
Lagrangian subspaces which are also solutions of (11.21) it will be useful to consider the 
simplest possible context for the Maslov index, linear systems on R 2 . To illustrate the role 
of A, a stability problem for a reaction-diffusion equation is used. It is a simplified version 
of the class of nonlinear parabolic PDEs studied in [Hj . 

Consider the nonlinear parabolic PDE 

dch d 2 6 



for the scalar- valued function (p(x,t). There is a basic steady solitary wave solution 

$(x) = | sech 2 , (5.2) 



9 



-e- 



-e- 



Fig. 1. Elements of the spectrum of Jzf in the complex plane. The big dots are the elements 
of the point spectrum Jzf and the dashed interval represents the essential spectrum. The origin 
A = is located at the intersection of the two axes. 



which satisfies 



0. Linearizing fl5.ll) about the basic state 4> and looking for 



solutions proportional to e leads to the spectral problem 



A0 , with Jz^> := 



dx 2 



(j) + 2(f)(x)(j). 



(5.3) 



The basic state (15.21) is said to be (spectrally) unstable if any part of the spectrum of Jzf 
is positive. The spectrum of Jzf can be explicitly constructed. It consists of a branch of 
essential spectra and a point spectrum 

a(Jz?) = a ess (jzf) U oy,(jSf) , 

with <r ess (Jzf) = {AgM : A < — 1} and cr p (Jzf) = j— § , , | j. The spectrum is illustrated 
in Figure [TJ 

The point spectrum can be verified by constructing the Evans function. First reformulate 
first-order system. Let 



u(x, A) 



<f){x,\) 



y(j) x (x, A)^ 



then 
with 



Ju x = B(x, A)u , ue 



A e 



(5.4) 



B(x,A) 



A + 1 - 3sech 2 (|x) 
1 



The eigenvalues of Aoo(A) are real and hyperbolic when A + 1 > 0. In this formulation 
the stable (u~) and unstable (u + ) subspaces are represented by 



u^A) =e ±7S 



i(^±7^), 



where s = ^x, 7 = 2\/\ + 1, 



^(s, A) = ±ao + aitanh(s) ± a2tanh 2 (s) + 03 tanh 3 (s) 
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and 



a = -^(4-7 2 )a 3 , a x = l(^ 2 - 3) a 3 , a 2 = a 3 , (5.5) 
15 5 



and a 3 is an arbitrary nonzero real number. The Evans function is then 

D(X) = u~(x,\) Au + (x,\) . 
Evaluating at x = 0, a straightforward calculation leads to 
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D(A) = A (4A + 3)(4A - 5) . 

The zeros of D(\) are the eigenvalues, confirming the point spectrum j — |, 0, +| |. 

For linear Hamiltonian systems on M 2 Lagrangian subspaces are just one-dimensional 
subspaces. The path of unstable subspaces u + (x, A) is used to define the Maslov index. 
The natural one-dimensional subspace to choose for the reference space is E^ C (A), 



E^(A) = span 



2 

-7, 



Then, assume simple intersections between E^(A) and u + (x, A) - which can be confirmed 
a posteriori for the example (15.41) - and assume that 

hm u + (x,A)nE^(A) = {0}. 

This latter assumption is equivalent to assuming that A is not an eigenvalue. The Maslov 
index for this case is 

Maslov(u + , ESJ = £ sign (Ju+ u + > vol , 

with xq the points at which u + (x, A) PI E^ is non-trivial, and the volume form can be 
taken to be vol = e± A This expression for the Maslov index is the one-dimensional 
version of (12. 7\i . 

The path of unstable subspaces is 



u + (x, A) = ^e 



( 2h~ 



(5.6) 
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The intersection form in this case is 



r(u+,E^,x ) = (Ju+ u+) 



vol 



X=XQ 



vol 



(«) 2 - A - 1 + 12sech 2 s)« 



vol . 



X = XQ 



However, at a point xq where u + intersects E A , u\ = —^uf and so 

T(u + , E^, x Q ) = 12 sech 2 |x «) 2 vol . 

Hence r(u + ,E 00 s ) > at each intersection, and the Maslov index is just the sum of the 
intersections. An intersection occurs when 



f Au + = 0, where £ s 



"7, 



Now 



( s A u + = 1u\ + . 
The factor e 7S is not important and so can be divided out, giving 

dh+ 



f A li- 



ds 



+ 2 7 /i + . 



This function has 0, 1,2 or 3 zeros depending on the value of A. Each zero corresponds 
to an intersection between the unstable subspace with E^(A). The function C, s A u + is 
illustrated in Figure [2] for the case A = —0.8 where £ s A u + has three zeros indicating 
three intersections. A summary of the Maslov index in each region is tabulated below. 



A 


— 1 < A < — | 


-f <A<0 


< A < | 




Maslov(u+,E^) 


3 


2 


1 






Note that the Maslov index jumps by one at each eigenvalue. Let Ao be any fixed real 
value of A such that Ao > — 1 and Ao is not an eigenvalue, then the value of the Maslov 
index equals the number of eigenvalues of ££ in the set A > Ao. 

This link between the Maslov index and the number of eigenvalues is not a coincidence 
specific to this operator but a quite general fact for systems such that <9aB(x, A) is a semi- 
definite matrix (see [15] for example). In the case where <9aB(x, A) is a not semi-definite, 
we give a counterexample in section [T71 
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Fig. 2. Plot of £ S (A) A u + (x, A) as a function of x for the case A = —0.8. 
5.1 The Maslov angle in ^(M 2 ) 

Another way to count intersections between the path u + and some reference plane is to 
use the Maslov angle (12.41) . In this case the angle k(x, A) is just the angle determined by 
a polar representation of u + 

e i K ( x ,x) . = ut(x,\)-iut(x,\) 
Ui(x, A) + i«2~(£, A) 

As x — > ±oo 

lim e iK{x ' x) = — - . 

The Maslov index is then the count of the number of times that k crosses some reference 
angle, such as the angle associated with the stable subspace. This version of the Maslov 
index is equivalent to the definition based on intersection index. 



6 The Maslov index for paths 

In this section we look at some of the properties of paths of Lagrangian subspaces that 
are also solutions of (11.21) . and bring in the exterior algebra representation. 

Here and throughout, let ei, . . . , be the canonical basis for M 4 . The symplectic form in 
standard form is then 

oj = ei A e 3 + e 2 A e 4 . (6.1) 
The symplectic form is related to the symplectic operator J in (11.31) by 

(Ju,v) = [w,uAv| 2l Vu,v6K 2 ". 
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Here and throughout (-, ■) is the standard inner product on M 4 and [■, is the induced 
inner product on A fc (R 4 )- We will sometimes write (-, -)d on IR d when the associated 
dimension is not clear. The equivalence between the induced inner product [•,•]& and 
(■, -)d when d = dim (^/\ k (M. 2n ^J is established in Appendix IDl The basic properties of 
symplectic exterior algebra can be found in [3"T] . 

The set of all Lagrangian subspaces of 1R 4 associated with the standard symplectic operator 
oj will be denoted by A(2). A(2) is a manifold of dimension 3. It is a submanifold of G 2 (IR 4 ), 
the Grassmannian of all 2— dimensional subspaces of R 4 [3"T] . 

A path of unstable subspaces will be represented by U + (x, A) € A 2 (R 4 )) which satisfies 
the equation (13.31) and f!3.4p . When A is not an eigenvalue then U + (x, A) — > E^ C) (A) as 
x — > ±oo and so in F(E u (x, A)), projective space based on E u (x, A), a loop of Lagrangian 
subspaces is obtained. The Maslov index is then taken to be the sum of the weighted 
intersections of U + (x, A) with E^ C (A). 

E^ C (A) is represented by the 2— form C _ (A) defined in ( 13. 6ft . An analytic basis can always 
be constructed for E^ (A) [8]. Denote this basis by 

E^(A) = span{^(A),e 2 s (A)}. 

Consider the 3— form 

U+(z,A)A(ai£(A) + a 2 g(A)). 

If for any fixed x and A, 

U+(x, A) A (c*i£(A) + a 2 &(A)) = a = 0, 

then we say that E^ C (A) is transverse to E"(x, A) := image(U + (x, A)) at that value of 
(a, A). 

We say that E^ C (A) and E u (x, A) have a simple intersection (or regular intersection) if the 
intersection is one dimensional; there exists a £ MP 1 such that 

U(x,A)A£ = with e = "^(A) + a 2 ^ 2 (A) . (6.2) 

We will assume that non-trivial intersections are regular. It is proved in [3] that non- 
regular intersections can be eliminated by perturbation. This property can also be proved 
using the homotopy equivalence property of the Maslov index, and a nice proof of this is 
given in §3.4 of [38J. 

Now, let uf(x, A), (x, A) be a basis for E u (x, A) such that 

U + (a;, A) = uf(x, A) A u^(x, A) . 
Then (E2D implies that there exists (/5i,/3 2 ) eR 2 \ {0} such that 

f = ft<(x, A) + /9 2 u+(x, A) := Z + (x, \)(3 . (6.3) 
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The Maslov index is a count of how many times the path XJ + (x, A) crosses the refer- 
ence space E^(A), weighted by the intersection form. Starting with (12.61) . and using the 
equation (II. 2p and the equivalence (I6.3p it is 

r(A,V,zo) = <B(z ,A)£,Ovol. (6.4) 

Note that the x— dependent path U + (x, A) drops out. For fixed A, once we have found a 
point xq corresponding to a regular intersection, the intersection form can be evaluated 
using eigenvectors spanning E^ C (A) and the known symmetric matrix B(xo, A). However 
there is the hidden calculation of determining a G MP 1 and an algorithm for this is 
developed in £ fT"""l 

Suppose A is not an eigenvalue. Then the Maslov index of the path U + relative to is 

Maslov(U + , E^) = £ sign (£, B(x , X)0 vol , (6.5) 

XQ 

where the sum is over all points x® of intersection in the interval — oo < xq < +oo. 
This expression will serve as a definition for the Maslov index of a solitary wave at A: 

Maslov(A) = Maslov(U + (-, A), E^(A)). (6.6) 

Proposition 1 Suppose A is not an eigenvalue. Then under the given hypotheses on Ql.ty) 
the Maslov index of a solitary wave is finite. 

Proof. Introduce a metric dist(- , •) on the manifold of 2— dimensional subspaces of 1R 4 . 
For example this can be the standard metric on the Grassmannian G2(M 4 ) [4*5] . 

Let AgK \ cr ess which is not an eigenvalue. Since E^(A) and E^ C) (A) are transverse, by a 
suitable scaling of the stable subspace E^(A) we can take 

dist((E^(A)),E^(A))>l, 

where (E^(A)) is the closed set of planes which are not transverse to (E^(A)). To simplify 
notation, the argument in dist(- , ■) should be interpreted as the representation of E s,u on 
the Lagrangian Grassmannian. 

When A is not an eigenvalue we have that E M (x, A) — > E^ Q (A) as x — > ±oo. Therefore, 
given any e > 0, there exists 5 > such that 

dist(E u (x, A), E^(A)) < e for \x\ > \ . 

o 

Now use the triangle inequality 

1 < dist((E^(A)), E^(A)) < dist((E s 0O (A)), E u (x, A)) + dist(E«(x, A), E^(A)) , 
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or 

dist(E tt (a;,A),(Ef (A))) > l-e for Id > - . 

o 

Hence there exists > such that for \x\ > E u (x, A) and E^ C (A) are transverse. 
Intersections are therefore limited to the finite interval — < x < x*. Since intersections 
are generically isolated their number is finite. ■ 



7 The Maslov index on A (K 4 ) 

The vector space A 2 (^ 4 ) is six-dimensional, and the orthonormal basis induced from the 
basis of M 4 is 



Ei = ei A e 2 , E 2 = ei A e 3 , E 3 = ei A e 4 , 

E 4 = e 2 A e 3 , E 5 = e 2 A e 4 , E 6 = e 3 A e 4 . 
Any U G A 2 (^ 4 ) ca n be represented in the form 

u = £^e,. 

i=i 

The Grassmannian Go(M. ) is a subset of A 2 (K 4 ) defined by 

= UAU = /ivol, h := UxU 6 - U 2 U 5 + U 3 U 4 . 

The Lagrangian-Grassmannian is the subset defined by 

= ui A U = I 2 vol , I 2 -=U 2 + U 5 . 

The Lagrangian-Grassmannian A(2) is the three dimensional submanifold of P (7\ 2 (I 
defined by I\ — J 2 = 0. 

Let V G A 2 (^ 4 ) be a fixed Lagrangian plane. Then 

A x (2) = { U G A 2 (M 4 ) n A(2) : U A V = } , 
is a codimension one submanifold of A(2) [3]. We have a sequence of manifolds 



(7-1] 



(7.2) 



(7.3) 



(7.4) 



Manifold 


A 2 (K 4 ) 


MP 5 


G 2 (M 4 ) 


A(2) 


A 4 (2) 


Dimension 


6 


5 


4 


3 


2 



In this table MP 5 represents P [A 2 (M. 4 ) 

Consider the class of linear Hamiltonian systems (11.21) with B(x, A) satisfying the asymp- 
totic properties fll.4p - fll.5D . 
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Proposition 2 A(2) is an invariant manifold of $3.3\) . 
Proof. 

£U A U = \] x A U + U A U x . 

= A( 2 'UAU + UAA< 2 >U 
= Trace (A) U A U 
= 

since Trace(A) = 0, also using the property [2] 

A (2) U A U + U A A (2) U = Trace(A)U A U . 
This proves that U A U is a constant along solutions. Similarly, 

£w A U = uj A U x 

= u> A A^U 

= u> A A^U + A^w A U - A^'w A U 
= Trace(A)w A U - A^uj A U 
= -A^u A U 

= 0, 

since ui is in the kernel of A^, a property which is proved in AppendixEl This proves that 
U A U and wAU are constant along solutions. Hence the special case UAU = o;AU = 
completes the proof. ■ 

Let V = span{£i, £2} be a fixed Lagrangian plane; that is £1 and £2 are linearly independent 
and (J£i,£2) = 0. The reference subspace V is represented by the form 

V = £iA&- 

The intersection between a path of Lagrangian subspaces TJ(x, A) and V can be described 
as follows. For each fixed A define 

gt{x) = {«6l 2 : U(x, A) A [axix + a 2 f 2 ) = o} . 

Then there are three cases 

• If = {0} then U and V are transverse at x = Xq. 

• If srf (x ) is one dimensional then U and V intersect in a one-dimensional subspace at 
x = xo (the case of regular crossing). 

• If g/(xo) = M 2 then U and V intersect in a two-dimensional subspace at x = xq (this 
case is sometimes referred to as an intersection between U and the vertex of V). 
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The Maslov index for a path of Lagrangian subspaces is given by (16. 5ft . However, in the 
case of A 2 (K 4 ) a new representation of the intersection form can be obtained. Suppose 
that a regular crossing occurs 

U n V = span{£} , 
at x = Xq, then the crossing form is 

r(U,V,x )=a,A£AA£, (7.5) 

To verify this formula, note that 

wAaA J~ c = (a, c) vol , for any a,ceK . 

Hence 

wA(AA( = wA(A J _1 Bf = (£, Bf ) vol , 
recovering the expression in (16 .4p . There is an interesting geometric interpretation of (17. 5|) . 
At a regular intersection the two-plane £ A A£ is not a Lagrangian plane. It is an element 
of G2(M 4 ) but not an element of A(2). Since A(2) is a codimension one submanifold of 
G2(M 4 ), the sign of u>A£AA£ determines which side of A(2) in G2(M 4 ) it lies. See Appendix 
[E]for further discussion of this case. 

Hence the Maslov index of a path U(x, A) relative to V is 

Maslov(U, V) = J2 si g n (^ A £ A A£) , 

where the sum is over all interior intersections. 



8 The Maslov angle on A 2 (^ 4 ) 

In M. 2 , the Maslov angle is just the angle associated with the polar representation of a 
vector in R 2 as shown in §5.11 For a Lagrangian frame of the form (12.11) the Maslov angle 
is defined as in (12. 3p and (12 .4p . In this section a new formula for the Maslov angle is 
given for the exterior algebra representation of a Lagrangian plane. Here, the result for 
4D phase space is given and in Part 2 [18], the general result for 2n— dimensional phase 
space is given. 

A Lagrangian frame can be partitioned into two 4x2 blocks as in (12 .ip and it can also 
be represented in terms of its columns: 




= [zi|z 2 ], with (Jzi,z 2 ) = 0. (8.1) 
Denote the exterior algebra representation of the Lagrangian plane by 

U = ZiAz 2 . (8.2) 
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Proposition 3 There exists a 2— form C, 

C = Ci + iC 2 , with d,C 2 eA 2 (K 4 ) 

such that 

det[X-iY]vo/= C AU. 



(8.3) 



It follows from this proposition that there exists a scalar complex- valued function K such 
that: 

C A U = AT(U)vol. (8.4) 

A formula for the Maslov angle n (the real number defined up to an addition by a multiple 
of 2tt by e lR = det((X - iY)(X + iY)^ 1 ).) is then immediate: 

Proposition 4 

e iK = AT(U)/F(U). 
It remains to prove Proposition [3j The proof is by explicit construction. Let 

c j = e j -i3e j , j = l,2. 

Then 

T 

/ {c l ,z 1 ) ( Cl ,Z 2 ) 
(c 2 ,zi) (c 2 ,z 2 ) 



X-iY 



V 



il 



X 
Y 



Ci C 2 Zi z 2 



(8.5) 



and so, using the induced inner produclO] on A 2 (K 4 ) (see Appendix [DJ 



det[X-iY]vol = det 
This gives a formula for K, 



(Ci,Zi) (ci,z 2 > 
(c 2 ,zi) (c 2 ,z 2 ) 



vol = [ci A c 2 , U] 2 vol 



K{U) = def Id Ac 2 ,U] 2 . 

It is not necessary to give an expression for C since in computation it is K that is needed. 
However, for completeness it is given. Let C be an 2— form satisfying 



Then 



d A c 2 A C = [d A c 2 , d A d| 2 vol . 



det[X-iY]vol = C AU 



(8.6) 



2 A real inner product is used throughout the paper. Complexification is used so rarely, a Her- 
mitian inner product is not necessary. One just needs to keep track of the complex conjugations. 
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The 2— form C is in fact the Hodge star of Ci A c 2 although the details of that character- 
ization are not needed. 

Now compute the formula in coordinates on A 2 (^ 4 )- On M 4 with the standard basis, 

Ci A c 2 = (ei - iJei) A (e 2 - iJe 2 ) 
= (ei - ie 3 ) A (e 2 - ie 4 ) 
= ei A e 2 — iei A e 4 + ie 2 A e 3 — e 3 A e 4 , 

and so, if U = Uj~Ej, with E 1; . . . , E 6 the standard basis on A 2 (^ 4 ), 

K(U) = [a A c 2 , U] 2 = U X - iU 3 + iU 4 - u 6 , 

and so the expression for the Maslov angle is 

= Ui-U*-iU 3 + iU* , 7) 

f/x-f/e + kVs-itV 1 ' J 

This expression is equivalent to the formula derived in equation (22) of [5]. 
The two-form C in (18.31) is computed to be 

C = -ei A e 2 + e 3 A e 4 + i(ei A e 4 — e 2 A e 3 ) . 

9 Further decomposition of the Maslov angle using the eigenvalues of a uni- 
tary matrix 

Let Z e R 4x2 be a Lagrangian frame on M 4 of the form (18.11) and U defined as in (18. 2|) . 
Then the matrix 

Q = (X-iY)(X + iY)-\ 
is a unitary and symmetric (but not Hermitian) matrix. 

The Maslov angle for a closed path (12.31) is defined using the determinant of Q. However, 
Q has 2 eigenvalues of unit modulus. Denote these eigenvalues by e 1Kj , j = 1,2 with kj 
real. Then 

e iK = e iKl e iK2 k = «i + k 2 (mod 2tt) . 

These eigenvalues are independent of the choice of Z e M 4x2 as a representative of a 
Lagrangian space: choosing another representation leads to similar matrix. 

These eigenvalues can also be used to give another formula for the sign of each intersection. 
Fix the reference angle to be (mod 2tt). There is a 1-dimensional intersection at Xq with 
a reference plane if and only if there exists e 1Kr ^ ' = 1 with r = 1 or r = 2. There is 
a 2-dimensional intersection at xq with the reference plane if and only if e 1Kr ^°- ) = 1 for 
r = 1 and r = 2. If e 1Kr ^ ^ 7^ 1 for r = 1 and r = 2 then the intersection is transverse. 
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When the intersection is regular, the sign of the intersection is given by: 



lim #{r G S\K r (x ) G (0, tt) + 2vrZ} - #{r G S\K r (x ) G (-7T, 0) + 2vrZ} . 



Thus, it is possible to determine the Maslov index, defined with intersections by simply 
tracking the crossings of the angles «j with 27rZ. 



9.1 The angles Kj in the exterior algebra representation 



The two angles k±, k 2 satisfy 



det(/xl - Q) = /i 2 - Trace(Q)/i + det(Q) = 



(9-1) 



with fi of unit modulus and fi r = e lftr , r = 1,2. Both the trace and determinant can be 
expressed in terms of the exterior algebra representation. 

For the determinant, as shown in 3HI 



det[X - iY] = K(U) := [ Cl A c 2 , U] 2 = U 1 - iU 3 + iU 4 - U e 



for A (M 4 ) 3 U = ELi UjEj. Hence 



det(Q) =K(U)/K(U). 



(9.2) 



It remains to express the Trace of Q in terms of the exterior algebra representation. 
Proposition 5 

2 



Trace (Q) 



K(U) 



(9.3) 



To prove this proposition, use ( 18. 5ft to relate the columns of Z to the X — Y decomposition 



X + iY 



(Ci,Zi) (Ci,Z 2 ) 

(c^,zi) (c^,z 2 ) 



Hence 



Q _l f ( Cl ' Zl ) ( C 1' Z 2) 

K \ (c 2 ,zi) (c 2 ,z 2 ) 



(c 2 ,z 2 ) -(ci,z 2 ) 

-(^Zi) (cT,Zi) 
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and so 



Trace(Q) = i((c 1)Z 



= i|det 



i) (c 2 ,z 2 ) - (ci,z 2 ) (c 2 ,zi) - (c 2 ,zi) (Ci,Z 2 ) + (c 2 ,z 2 ) (Ci,Zi)) , 
det 



(Ci,Zi) (ci,z 2 ) 
(ci,zi) (c£, z 2 ) 



(Ci,Zi) (Ci,Z 2 ) 

(c^,zi) (ci,z 2 ) 



= i ([ci A c 2 , Zi A z 2 ] 2 + [ci A c 2 , Zi A z 2 ] 2 } 

= i ([ci A ci + cT A c 2 , zi A z 2 ] 2 ) 
= | ([ex A e 2 + e 3 A e 4 , z x A z 2 ] 2 ) 
= |(£/ 1 + f/ 6 ), 



using 

Re (ci A ci) = Re ((ei - ie 3 ) A (e 2 + ie 4 )) = ei A e 2 + e 3 A e 4 , 

proving (19.3j) . 

Given a path U G A 2 (^ 4 ) the eigenvalues of Q can be computed by substituting (19.21) 
and (E3} in (EH) leading to 



1«1,2 



<7i + ?7 6 ± v^lWe + 2LW 4 - Uj - Uj 
Ux + iU 3 - iU 4 - U 6 



Using the properties of a Lagrangian plane, U 2 + U 5 = and Z7i?7 6 — U 2 U 5 + L^L^ 
this formula reduces to 



11 ■— o 1K l,2 

A*i,2 •— e 



[7 1 + ^ 6 ±i v /4^ 5 2 + (^3 + ^4) 2 

c / 1+ i[/ 3 _i[/ 4 _ t / 6 



The formulaGD for the decomposition of the Maslov angles in dimension 3 is given in Part 

2 nsi. 



3 In the most general case, the eigenvalues e 1Kj are the roots of the following polynomial: 

P(A) = det((X - 1Y) - /x(X + iY)) , liGS 1 . 

The coefficients of P are antisymmetric multi-linear functions of Z. As a consequence, they can 
always be expressed as a linear combination of the minors of Z and hence, of the components 
of U. Therefore, it is possible to compute the eigenvalues e lKj from U. 
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10 A dependence of the Maslov index 



Let B(x, A) be as defined in ( 11. 2ft with the asymptotic property ( 11. 5p . When the A— dependence 
of B(x, A) takes a simple form one can say more about the A— dependence of the Maslov 
index. For example, in §|H]the matrix B(x, A) in (15. 4p satisfies 



| B (,,A) 



1 




that is, it is positive semi-definite. In the example in 
with A(x, A) defined in ( I11.3P has the property 



the matrix B(x, A) = JA(x, A) 



fjB(*,A) 



-10 
0-100 





which is negative semi-definite. 

When d\B(x, A) is semi-definite, the Maslov index is a monotone function of A. For the 
case of gradient systems (as in §0 and ^TTT) this property is proved in Lemmas 3.3 and 3.7 
of Bose & Jones [6]. Related results are proved by Arnold [I] and generalizations of 
these results are proved by Chardard [15]. These results are summarized in 

Lemma 6 Assume that: 

• B(x, A) is a smooth function with respect to x and analytic with respect to A. 

• There exists B^A), 7 > and F > such that\/x,X ||B(x,A) -B oc (A)|| < Fe~^ xl . 

• The open set X = R — a ess of real numbers is not empty. 

• d\B(x, A) is semi-definite symmetric matrix. 

If [Ai, A2] fl a ess = and Ai, A2 ^ cr, then Maslov(A2) — Maslov(Ai) is equal to the number 
of eigenvalues with multiplicity in [Ai, A2]. 

The first three assumptions are the usual hypotheses made to prove the analyticity of the 
Evans function and the theorems linking eigenvalues and the zeros of the Evans function. 
Using this lemma, it is possible to define the following invariant for the homoclinic orbit: 

Definition 7 MasW lomoclinic is defined as lim A _ + Maslov{\). 

If we make the hypotheses of Lemma El Maslov homoclimc is only dependent on A(x, 0), and 
hence on the linearization of the ODE satisfied by the homoclinic orbit. 



23 



In fact, it is possible to define Maslov homoclimc without any referenced to a parameter A 
(see [5115115] ). In these references, the quantity Maslov homoclinic is defined like a Maslov 
index for the path Z(., 0) with respect to E s (0). Unfortunately, Z(., 0) does not necessarily 
admit a right end point, and when it has one, the intersection with E s (0) is not trivial 
and some extra work is needed. 

Anyway, these geometrical constructions are not convenient for numerical computations 
and we will not use them in the sequel. 

10.1 The Maslov index for large values of \X\ 

When A — > +oo (or A — > — oo if the essential spectrum extends to minus infinity) we 
expect the Maslov index to converge to some finite value. This property is similar to the 
property of the Evans function for large A. The hypotheses are based on the analogous 
result of jlS], adapted to the setting of the Maslov index. 

Hypothesis 8 Suppose 

• that there exists A el such that a ess is empty for all A > A ; 

• For large enough X, A^A) has no purely imaginary eigenvalues; 

• Let K(A) be a symplectic 4x4 matrix depending analytically on A whose first 2 columns 
are a basis for E^(A) and whose last 2 columns are a basis for E^ C (A) Define 

F(x, A) = K~ 1 (A)(A(x, A) - A oc (A))K(A) . 

• Suppose that, for large enough A: 

/ R \F(x, X)\dx is bounded, uniformly in X 

I\x\>x |F(x, X)\dx tends to when x — > oo, uniformly in X 

f R |F( 2 )(x, A)ei|dx tends to . 

Remark. If there exists Ao such that a ess is empty for all A < Ao the above hypotheses 
can be modified accordingly. 

Proposition 9 Assume that hypothesis^ is met by A(x, A), then 

lim D(X) = 1 , and lim Maslov(-; A) = . 

A— >+oo A— >+oo 



For example, in the case of a transverse orbit (i.e. dim(E s (x,0) n E u (x,0) = 1), we have: 

Maslov homoclinic = lim^ +00 lim e _ + £>0 max( Maslov(Z(] - oo,x - e], 0), E s (x, 0)), 

Maslov(Z(] - oo, x + e], 0), E s (x, 0))). 



24 



The proof is obtained by following the argument in Proposition 1.17 in Pego & Wein- 
stein |3S] and the Appendix of Bridges & Derks [S]. Let JC(\) be the matrix whose 
entries are the 2x2 minors of K(A). In terms of bialternate product, /C(A) = K(A) ©K(A). 
One proves that )C(\)~ 1 Y(., A) converges, uniformly in x, to the constant vector ei when 
A — > — oo. Then, for large enough A, /C(A) _1 F(., A) has a null Maslov index and so does 
Y(.,X). See pj)] for a detailed proof. 

10.2 Defining a Maslov index at A = when the basic state is approximated by hyperbolic 
periodic solutions 

In this paragraph, we assume that B(x, 0) = D 2 H(<f)(x)) where H : IR 4 — > R is the 
Hamiltonian function and is the basic homoclinic solution of the autonomous system 

Ju x = V#(u). (10.1) 

This hypothesis is natural in the solitary wave context and is indeed satisfied by all 
examples presented in this article. 

Suppose that d\B is semi-definite near A = and that the dimension of the space of square 
integrable solutions of Ju x = B(x, 0)u is one. Furthermore, suppose is approximated 
by hyperbolic ^-periodic solutions fa, for which the Maslov index is well-defined [TofTl] . 

If the Maslov index at A = is defined as the limit of the Maslov index of the periodic 
solutions fa, when k — *■ 0, then, as shown in [13] under natural hypotheses, the Maslov 
index at is the value of limit of the Maslov indices of the periodic orbits when A is close 
to and it has the sign of f'(k) near k — 0, where f(k) := if 



11 A coupled reaction-diffusion equation with explicit Maslov index 

Consider the system of react ion- diffusion equations 

du 
~di 
dv 
dt 

where c, the coupling constant, is a non-zero real parameter, restricted to the values 
c > —2. When c > —2, it is straightforward to show that the trivial solution u = v = is 
stable in the time dependent problem, and the trivial solution of the steady equation is 
hyperbolic. 

This system has the exact steady solitary-wave solution 

u = v := u(x) = sech 2 (a;) . 



d 2 u 
dx 2 

d 2 v 
dx A 



Au + 6m 2 — c(u — v) 



— 4v + 6v + c(u — v) , 



11.1^ 
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Linearizing (lll.ip about the basic state u and taking perturbations of the form 

e xt (u(x,X),v(x,X)), 
leads to the coupled ODE eigenvalue problem 

u xx = (A + 4 + c — 12u(x)) u — cv 

v xx = — cm+(A + 4 + c— 12u{x)) v . 

This eigenvalue problem can be written in the standard form 

u x = A(x, A)u, u e M 4 , 

with u = (u, v, u x , v x ) and 

( 1 (A 
1 

f(x, A) -c 
v -c /(x,A)0 0y 

The system fl 11.31) is Hamiltonian: JA is symmetric. 



;ii.2) 



;n.3) 



A(x,X) 



with f(x, A) = A + 4 + c-12 sech 2 (x) 



The spectral problem All .3j) can be solved explicitly. Write the second-order problem in 
the form, 



u 



V 



1 

1 



;ii.4) 



Let 



u = u — V 



V = U + V . 

Then substitution into (jll.4p leads to the decoupled system 

u xx + 12 sech 2 (x) u = (A + A)u 
v xx + 12 sech 2 0) v = (A + 4 + 2c)v . 



;il.5) 



These two systems have explicit solutions (cf. Appendix [F]), and using these results one 
finds that there are exactly six eigenvalues for the spectral problem (111.21) : 

Ai = — 3 — 2c , A2 = — 3 , A3 = —2c , 
A4 = , A5 = 5 — 2c , Ag = 5 . 
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Fig. 3. Plots of the spectrum in the A-plane for c = —2 (starting at the top), c = —1, c = 0, 
c = 1 and c = 3. For each plot, the intersection of the axes is the origin A = 0. 

The essential spectrum is 

a ess = {A G K : A < —4} U{AGl : A < -4 - 2c} . 



When c = then there are three double eigenvalues: 

Ai = A2 = —3 , A3 = A4 = , A5 = A6 = 5 . 



For c nonzero and small the eigenvalues Ai, A3 and A5 are perturbed to the left (when 
c > 0) and to the right (when c < 0). Hence positive coupling is stabilizing and negative 
coupling is destabilizing. 

When c = —2 (the lower bound on c for stability of the zero state and hyperbolicity) 
there are four positive eigenvalues, which is the maximum number of positive eigenvalues. 

At c = — I one of the positive eigenvalues passes through zero, leaving 3 positive eigenval- 
ues. Then when c = another eigenvalue arrives at zero leaving two positive eigenvalues. 
Then when c = | a third eigenvalue passes through zero. For all c > | there is one positive 
eigenvalue, one zero eigenvalue and four negative eigenvalues. The configuration of the 
eigenvalues as a function of c is shown in Figure El 

When the system decouples into two subsystems, the Maslov index is the sum of the 
Maslov indices of the subsystems 



Maslov 2D © Maslov 2 - = Maslov 4D 



;il.6) 



If we choose a Lagrangian frame of the form 



x 2 

Yx 



to represent the unstable 
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space, we have indeed det(X + iY) = det(Xi + 1Y1) det(X 2 + iY 2 ). From this property, 
(111.61) follows easily. 

This property is obvious in the present example, since the system decouples. The Maslov 
index for the first 2D system in (111.51) is tabulated below. 



A 


-4 < A < -3 


-3 < A < 


< A < 5 


A > 5 


Maslov^ 


3 


2 


1 






The Maslov index for the second 2D system in (111.51) is tabulated here. 



A 


-4 - 2c < A < -3 - 2c 


-3 - 2c < A < -2c 


-2c < A < 5 - 2c 


A > 5 - 2c 


Maslovif 


3 


2 


1 






The Maslov index for the full 4D system for any A is then obtained by fixing c and then 
applying the sum formula (111.61) . For example, fix c = — 1 and compute Maslov homoclimc , 



Maslov 



homoclinic 



c=-l 



lim 



Maslovf © Maslov 



2D 

■"> 



3. 



c=-l 



A summary of the Maslov index of the homoclinic orbit as a function of c is given in the 
following table. 



A 


-2<c<-§ 


-§ < c < 


< c < | 


of 


Maslov homoclinic 


4 


3 


2 


1 



12 Numerical implementation — approximation by periodic orbits 

When the solitary wave is approximated by a hyperbolic periodic orbit, the Maslov index 
of the periodic orbit is computed using the Maslov angle (I2.4p -( l2~5j) in the exterior algebra 
representation (e.g. equation (18.71) ). An algorithm for this case has been proposed in |16j . 
and a proof of convergence of the Maslov index in the limit as the periodic solution 
converges to the solitary wave is given in [H] . A first possibility would be to compute the 
Maslov index for a sufficiently good periodic approximation of the solitary wave and to 
compute the its Maslov index. 

However, it is possible to avoid the computation of periodic approximants and to adapt 
the algorithm presented in [16] to the homoclinic case. Fix A and a solitary wave. The 
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steps in this new algorithm are as follows. 

(1) Choose a large enough interval [—L,L]. 

(2) Compute the eigenvalue with largest real part of A^'(A), denoted by cr+(A), and its 
associated eigenvector C + (A)- 

(3) Integrate equation 

U+ = [A( 2 )(x,A)-a + (A)I]U + , (12.1) 

on [-L, L], taking C + (A) as initial condition at x = —L, using any standard numeri- 
cal integration scheme. The justification for the arbitrariness in choice of numerical 
scheme is given in Appendix 

(4) U + (L, A) and U + (— L, A) are nearly collinear, and an approximation to the Evans 
function is determined from their constant of proportionality 

U + (L, A) = D(A)U + (-L, A) + Error , 

where the error is generally of the order of machine precision. 

(5) Compute e 1 *^ using equation (18.71) or analogous representation. 

(6) Compute a lift of k(x) and choose the stepsize Ax so that \k(x + Ax) — k(x)\ < ir. 

(7) Compute the Maslov index using (12.51) . 

There are a number of sources of error in the algorithm. There is an approximation error 
due to the fact that the solitary wave is approximated by a periodic orbit. Two parameters 
have to be chosen: L and the step size Ax. The choice of step size is a familiar source 
of error. The consistency error of the numerical integration scheme will be of the form 
C Ax p , for some natural number p, at each step, where C is a constant depending on the 
derivatives of A. The choice of numerical scheme will also impose some stability condition. 

Since the Maslov index is an integer, the proposed scheme will give the Maslov index if 
the relative error on U + (-, A) is small enough. However, if sup R ||U + (-, A)|| is very small 
(for example when the Evans function is small), the relative error may be too big and 
lead to a miscomputed Maslov index. This is the case when A is an eigenvalue or near an 
eigenvalue, since there is an integer-valued jump in the Maslov index at eigenvalues. 



13 Numerical implementation — intersection index algorithm 

The numerical algorithm based on the intersection index is similar to the algorithm in 
£fT2l except that the computation of the angle k(x) is replaced by the computation of the 
angles K\ and n 2 . The previous algorithm can be modified as follows: 

(1) Choose a large interval — L < x < L. Initialize Maslov to 0. 

(2) Construct a symplectic matrix K(A) such that 



K(A) 
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represent the stable and unstable spaces of Acq (A). K(A) defines a symplectic change 
of coordinates in which the coordinates of stable and unstable spaces in the exterior 
algebra are respectively 



U 







and Vn 



o 

V 1 / 



Let /C(A) be the matrix whose entries are the 2x2 minors of K(A). In terms of 
bialternate product, /C(A) = K(A) K(A). 

(3) Compute the eigenvalue with largest real part, o + (A), of A^(A). 

(4) Define U + (x, A) = JC(\)~ l XJ + (x, A) and integrate the equation for U + , 



d_ 

dx 



U+ = [/C(A) _1 A^(x, A)/C(A) - a+(A)I])U+ 



(13-1) 



-L, using any standard 



on [—L,L], taking Uo as initial condition for U + at x - 
numerical integration scheme. 

(5) Compute the angles {k\ and K2 corresponding to U + over [—L,L]. If an angle Ki 
crosses 2-7rZ between x and x + Ax, update the value of the Maslov index to: 

Maslov i — > Maslov + sign (/^(x + Ax) — Kj(x)) . 

(6) Return U + (L, A) A Uo as an approximation to the Evans function. 

(7) At x = +L, return the value of the Maslov index. 



14 The Maslov index of solitary wave solutions of KdV5 

In this section we study the Maslov index as a function of A for the ODE eigenvalue 
problem 

4>xxxx - P<Pxx + a(x)<f) = \<p , (14.1) 
where 0(x, A) is scalar valued, P is a real parameter and a(x) is a localized function 
which satisfies — » ±oo, with exponential decay of a(x) at infinity. For 
definiteness it is assumed that > 0. The ODE (114.11) can be put in the form (11.21) with 



u 



\ 



and B(x, A) 



a(x) - A 

-10 

1 

1 P 



;i4.2) 
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The spectrum of the system at infinity Aoo(A) = J 1 B oc (A) has the characteristic poly- 
nomial 

det[A oc (A) - fil] = - P/j, 2 + - A , (14.3) 
With > and A = the four roots are hyperbolic for all P such that 

P + 2^d~> 0, 

which is assumed to be satisfied henceforth. When A 7^ the essential spectrum will form 
the boundary of the hyperbolic region. The essential spectrum is 

a ess = {A E R : A = floo + Ps 2 + s 4 , s G R} . 

When 

A<A cdgc = a 00 -^P(P-|P|), 

o 

the spectrum of A 0O (A) is hyperbolic. Hence, all the hypotheses for the existence of the 
Evans function and the Maslov index are satisfied. We will apply this theory to determine 
the Maslov index of a class of homoclinic orbits. 

The eigenvalue problem (jl.2p appears in the linearization about a solitary wave solution 
of the fifth-order Korteweg-de Vries equation (KdV5). KdV5 appears for example as a 
model equation in plasma physics, and in the study of capillary-gravity water waves 
[28U3nif2drT8lf7Tf23Tf9] . 

To see the role of fll.2f) in the linearization of KdV5, consider the following form of the 
fifth-order KdV equation relative to a moving frame, moving at speed c, 

dt ° dx 9i v ) dx 3 dx 5 ^ ~ ^ ^ 

A further scaling can be introduced so that c = 1, but including c is useful for comparing 
with results in the literature on KdV5. Effectively, q is a third parameter, but its value is 
restricted to natural numbers. 

Steady solutions of ( 114.4[) . that decay to zero as x — > ±00 satisfy the 4 th — order ordinary 
differential equation 

- P^ XX + c0 - = 0. (14.5) 

The system (114.41) linearized about a solitary wave <j)(x) solution of (j!4.5p takes the form 

with 

^(f> ■= (j>xxxx -P(t)xx + c<P-(q+ l)$(x) q <P . (14.6) 
In this c. There are two spectral problems: 

= A0 and L0 = A0, L0:=^-^f. (14.7) 

dx 
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The operator J?f is self-adjoint (in a suitably-chosen Hilbert space) and so A G 1 whereas 
L is not self-adjoint and A - which is the stability exponent - can in general be complex. 
The relationship between these two eigenvalue problems is discussed in §151 First the 
Maslov index of the spectral problem Jzf0 = A0, which can be put in the form ( 11. 2ft . is 
studied. 

The ODE ( 114.51) has been extensively studied and many solitary wave solutions have 
been found; a classification is given in [TT]. There are some special cases where explicit 
solitary wave solutions can be constructed. An example is the explicit solution <f)(x) = 
||sech 4 (rr7g) which exists when q — 1, c = 1 and P — ^. However, the interesting 
solutions of ( 114. 5ft need to be computed numerically. They can be computed using a 
spectral method (approximate the solitary wave by a periodic function of large wavelength 
and then use Fourier series to represent it), or in the case of symmetric solitary waves a 
shooting algorithm can be used. We used both methods to compute solitary waves. An 
example of the family of one-mode solitary waves as a function of P, computed using a 
spectral method, is shown in Figure HI 

Although these solitary waves are solutions of the model ODE, they are representative 
of solutions of the full water-wave problem. Dias, Menasce & Vanden-Broeck [M] 
have found large-amplitude branches of these solutions in the full water-wave equations. 

Symmetric solutions are computed numerically using a shooting method: the starting 
point is an element of the tangent space of the unstable manifold and the ending point 
is a symmetric point. The time step typically used is and the integrator chosen was 
the fourth-order Runge-Kutta method. For this integrator and this time step, the error 
on the Evans function at (which is then equal to 0) was of order 10 -12 . This is quite 
close to the minimum error reachable by using double precision. 

The ODE (114.51) can be characterized as a Hamiltonian system on M 4 and the Hamiltonian 
in the original coordinates is 

E{<f>) = ¥1* + \ p <i>l - W + -^# +2 - M*** ■ ( 14 - 8 ) 

q ~\~ z 

The function E{<p) is constant along solutions of (I14.5f) (i.e. 4^ = 0). Physically, for 
equations like KdV5, this quantity is associated with momentum flux. For simplicity, 
we will just refer to it as energy. The energy of the periodic approximants gives some 
information about the nature of the limiting homoclinic orbit. 

The energy can be plotted as a function of wavenumber k along a branch of periodic 
solutions as the wavelength tends to infinity (k — * 0) as a function of q and P. First the 
case P = y , c = 1 and q = 1 is considered and it is shown in Figure El 

In this case the energy is a monotone function of wavenumber and the convergence k — > 
is rapid. Keeping q — 1 and c = 1 but decreasing P to P = — 1 begins to show oscillations 
indicative of a Shilnikov-type bifurcation as shown in Figure EJ Decreasing P further to 
P = —1.9 shows more dramatically the Shilnikov-type oscillations, as shown in Figure 
[3 In Figure [7] a sequence of bifurcations occurs along the branch. Each point on the 



32 



Fig. 4. Numerically computed solitary waves for the KdV 5 equation for the case q = 1, c = 1 
(ooo = 1) and — 2 < P < 3. For each value of P, the value of the corresponding unimodal 
homoclinic solution (f> is plotted function of x. 



energy- wavenumber diagram where E'(k) = corresponds to a saddle-centre bifurcation 
of Floquet multipliers. There are always two Floquet multipliers at +1 due to the fact that 
(114.51) is autonomous. When E'(k) = an additional pair of Floquet multipliers coalesces 
at +1. Each one of these saddle-centre bifurcations of the branch of periodic orbits leads 
to a secondary homoclinic bifurcation [TU]. So, in addition to the limiting homoclinic orbit 
that we are principally interested in, there is a countable number of other orbits generated 
along the branch, which are homoclinic to the branch of periodic orbits. Although there 
is an infinite number of bifurcations along the branch our numerical results show that the 
Maslov index of the limiting homoclinic orbit is finite. 
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0.1 0.2 0.3 0.4 0.5 0.6 

Wave number k 



Fig. 5. Energy of the ^-periodic solutions as function of k for q = 1, c = 1 and P = ^r. 



0.15 - 
0.1 - 
0.05 - 




Wave number k 

Fig. 6. Energy of the ^-periodic solutions as a function of k for q = 1, c = 1 and P = — 1. 
14-1 Computing the Maslov index as a function of A 

First consider the case P — q = 1 and c = 1 where the unimodal solitary wave is given 
explicitly. The lifts k(x) of the Maslov angle for this system are plotted as a function 
of x in Figure \E\ for various values of A. In Figure [9J the corresponding Maslov indices 
have been plotted as a function of A. The Evans function shows that J?f has exactly three 
eigenvalues in this case. Denote these eigenvalues by 

Ai < A 2 = < A 3 . 
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0.05 0.1 0.15 0.2 0.25 

Wave number k 



Fig. 7. Energy of the ^-periodic solutions as a function of k for q = 1, c = 1 and P = —1.9, for 
low values of k. 




Fig. 8. k as a function of x for the following values of A:— 2, —1.9, . . . ,0.9. k is A-growing. The 
parameter values are P = c = 1 and 3 = 1. 

The qualitative behaviour of the Maslov index in this case is similar to the example on M 2 
in £j5j The values of the Maslov are shown in the table below. The Maslov index in this case 
is computed using the Maslov angle, and this Maslov index is denoted by Maslov(/t, A). 
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Maslov index 
3.5 F ' 

3 - 
2.5 - 

2 - 
1.5 - 

1 - 
0.5 - 

i - 




Fig. 9. Evans function and Maslov index as a function of A for the explicit unimodal solitary 
wave solution when P = c = 1 and q = 1. In this case, the Maslov index is A-growing. 



A 


A < Ai 


Ai < A < A 2 


A 2 < A < A 3 


A> A 3 


Maslov(K, A) 





1 


2 


3 



Note that the Maslov index in each region predicts the number of eigenvalues of «5f in 
each A interval. 



A region 


A < Ai 


A < A 2 


A < A 3 


A < A edge 


# Eigs(Jgf) 





1 


2 


3 



It is immediate from this table that 



Maslov homoclinic = lim Maslovf/t, A) = 2 . 

The operator J*f has exactly one negative eigenvalue in this case. Our calculations indicate 
that this is the case for all the unimodal homoclinic orbits. It is easy to show analytically 
that the Maslov index of a unimodal homoclinic orbit is greater than or equal to 2. An 
elementary proof in given in Appendix [HJ This result has implications for stability of the 
solitary waves as solutions of KdV5 and it is discussed in £ fl~5l 

To test how accurately the Lagrangian Grassmannian is preserved by the numerical 
scheme, the values of 

Ji = U A U and / 2 = wAU, (14.9) 
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Fig. 10. The decomposability of a 2-form is equivalent to l\ = where I\ is defined in (114. 9D , 
The upper figure shows the logarithm of the difference between a(x) and its limit aoo, and the 
lower figure shows the logarithm of the value of I\ when q = 1, P = Ax = .01 and A = — 10 
The fast and irregular oscillations are associated with round-off errors. 



are computed as a function of x. In these calculations the standard explicit fourth-order 
Runge-Kutta algorithm is used. The value of I\ is shown in Figure [TU] and shows that the 
error is of order of the machine accuracy, except for a small region around zero, but the 
error there is still exceptionally small. Concerning I2, it is in fact exactly preserved, even 
numerically: if I2 is below machine accuracy, then it remains so. 

Using Proposition [9] and the proof in Appendix [H it follows that the KdV5 system satisfies 
hypothesis El and therefore -D(A) — > 1 and the Maslov index tends to as A — ► —00. 



15 Spectrum of Jz? and the stability of solitary waves for KdV5 

One of the intriguing properties of the Maslov index is the connection between the number 
of eigenvalues in subsets of the A space, the Maslov index and the stability of solitary 
waves. We know that the Maslov index counts the eigenvalues of =Sf and it turns out that 
this number plays a role in previous stability results. 

For KdV5 the connection between stability and the Maslov index is not obvious. For 
unimodal solitary-wave solutions of KdV5, Kodama & Pelinovsky [SU] have studied 
this connection and they show the following result. Suppose the following integral exists 

N(c,P) = [ c, P) 2 dx , 
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and is a differentiable function of c, and define 



if M(i jP ) <0 

1 if f (1,P)>0 ' 

The functional N(c, P) is sometimes called the momentum of the solitary wave. If q — 1, 
it is proved by Lewandosky [33] that r = 1 for all admissable P. 

In [30] it is argued (see proposition 3.8 there) that a unimodal solitary wave is stable if 
r = +1 and #Jz?~ = 1, where is the number of negative eigenvalues of . 

This observation is consistent with the theory of [7J where a instability results for a class 
of unimodal solitary waves were obtained. 



15.1 Stability of two-pulse solitary waves 

More refined results on stability of two-pulse solitary waves were obtained by Chugunova 
& Pelinovsky [20]- Suppose that L has only simple eigenvalues except a double eigen- 
value at and suppose ^(1,P) ^ 0. In [2D] it is proved for symmetric solitary wave 
solutions that 

N uns t = — r — N imag , 

where N unst is the number of eigenvalues with strictly positive real part of L and N~ mag is 
the number of pure imaginary eigenvalues of L with negative Krein signature. Using the 
formula #j£f~ = Maslov homoclinic - 1 gives 

N unst = Maslov homoclinic - 1 - r - N~ mag . 

It is immediate from this formula that if N~ mag = and Maslov homoclimc = 2 then the basic 
state is stable if ^ > and unstable if ^ < 0. 

dc dc 

Using the classification of Buffoni, Champneys & Toland [II], a two-pulse solitary 
wave has the classification 2(£) where £ is a natural number. In [17] . it is found numerically 
that two-pulse solutions have the following formula for the Maslov index 

■ 3 if £ is even 

Maslov homoclimc = { . (15.1) 

4 if £ is odd 

We can make several observations using this formula for the Maslov index of two-pulse 
homoclinic orbits. Suppose r = +1, then 

■ 1 — N~ rinn if £ is even 

Nunst = Maslov — 2 — TV, 



imag 



2-N- mag iff is odd 
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We have the immediate observation that a necessary condition for a 2— pulse homoclinic 
orbit to be stable is N^ mag > 0. From the parity of N^ mag) we have N unst = 1 and N~ mag = 
when I is even and the solitary wave is unstable. When I is odd, the parity of is not 

sufficient to determine the stability of the solitary wave. 

Buryak & Champneys [12] used a completely different method to study stability and 
they found that 2— pulse solitary waves are stable if £ is odd and unstable if £ is even, 
assuming that r — 1. This is consistent with the value we found for the Maslov index. 
When I is odd, the stability of the solitary wave is equivalent to N^ mag = 2. 

The Maslov index does not give any information about the purely imaginary eigenvalues, 
and so to determine their number, a calculation of the spectral problem is necessary. 
Some results on this are reported by Chardard [T3]. There it is found that when I is 
odd, there are eigenvalues with non-zero imaginary part, but they appear to have very 
small real parts. Further results are necessary to be certain about the spectral stability 
of 2— pulse solitary waves when I is odd. 

One way to check whether the real part of a complex eigenvalue is nonzero is to use the 
formula (for the case q = 1) 

r+oo ^ 

Re(X) = ----, <j) x \u(x)\ 2 dx . (15.2) 

\\U\\ Z J-oo 

Here A is the complex eigenvalue associated with the stability exponent and u is the 
associated eigenmode: 

^ (uxxxx - PUxx + CU-2(j)Uj = \u. (15.3) 

The formula (I15.2[) is derived by multiplying (115. 3p by the complex conjugate of m(-) and 
integrating over R: 

A||w|| 2 = / \u\ 2 cf> x + [ 4>u x ll= [ \u\ 2 (p x + [ ^\(u x u-uu^)= [ \u\ 2 <p x + i / 4>lm(u x u) 
Jr Jr Jr Jr I JR JR 

In spite of the simplicity of the formula (115.21) there is not much that one can say in 
general. If 4>(x) is an even function then <j) x is an odd function. Then it is immediate that 
|it(x)| 2 even implies that Re(A) = 0. However, this is a highly special case. 



16 A model PDE for long- wave short-wave resonance 

In this section the Maslov index is computed for a class of solitary waves which arise in a 
model PDE for long- wave short-wave (LW-SW) resonance (cf. Kawahara et al. [29], 
Ma [35], Benilov & Burtsev [5], Latifi & Leon [32]). The LW-SW equations are a 
coupled system with one equation of nonlinear Schrodinger type and the other of KdV 
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type. A typical form is 

E t = i(E xx + pE- vE) 

K H 1 (16.1) 

pt = d x (p xx -cp + 3p 2 + \E\ 2 ) , 

were p(x, t) is real valued, E(x, t) is complex valued, and c, v are considered to be positive 
real parameters. In real coordinates, E = u + iv and p = w, the above equations can be 
written: 

u t = —v xx — vw + vv 

Vt — u xx + uw — uu (16.2) 

Wt = Wxxx — cw x + 6ww x + 1uu x + 2v v x . 

This system can be expressed as a Hamiltonian system in the time direction. However, we 
will not emphasize this property since it is the spatial Hamiltonian structure that is asso- 
ciated with the Maslov index (see Appendix [J] for the temporal Hamiltonian structure). 
Solitary waves satisfy the steady equations 

— 2u xx — 2uw + 2vu = 

-2v xx - 2vw + 2uv = (16.3) 
—w xx + cw — 3w 2 — u 2 — v 2 = constant , 

where the signs and coefficients are modified to ensure that they are the Euler-Lagrange 
equation associated with the Hamiltonian function H(Z) in Appendix [J] Exact solutions 
of this problem are known [33]; for example, 

u(x) = A sech(y/u x) , v(x) = and w(x) = 2z/sech 2 (y / ^£) , (16.4) 

with constant = and A 2 = 2v[c — 4v), and the existence condition c — 4u > 0. 

To study the Maslov index of these solutions, linearize the steady equations about the 
basic solitary wave and introduce a spectral parameter: LZ = XZ with L = D 2 H(Z). 
Written out, this equation is 

—2u xx — 2wu — 2uw + 2vu = Xu 

-2v xx - 2wv - 2vw + 2vv = Xv (16.5) 
—w xx + cw — Qww — 2uu — 2vv = Xw 

When v — this system decouples into a second order equation for v, and a fourth order 
coupled system for u,w, 

—2u rr — 2wu — 2uw + 2vu = Xu 

(16.6) 

—w xx + cw — 6ww — 2uu = Xw 
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The decoupled equation for v is then 



2v„ 



2wv + 2vv = Xv 



(16.7) 



This latter system can be analyzed completely and the result in given in Appendix iKl 

The fourth-order system for u, w (116. 6ft can be written as a standard Hamiltonian ODE 
in the form (II. 2\i with n = 2 by taking 



u(x, A) 



/ u \ 



w 

2u,,. 



B(x,A) 



V 



A - 2v + 2w(x) 
2u 





2u(x) 







A - c + 6w(x) 

ho 





2 

1 



\w x J 

The essential spectrum for this equation is 

a ess = {AelR : \ >2v and A > c} . 

Adding the condition that c > 4u, the system at infinity is hyperbolic for all A e M such 
that A < 2v. 

The Maslov index of the 4x4 subsystem is computed for the case c = 1 and v — 0.2 
and the results are shown, along with the associated Evans function, in Figure [TT] and 
tabulated in the table below, where Ai < A2 = < A3 are the three roots of the Evans 
function. 



A 


A < Ai 


Ai < A < A 2 


A 2 < A < A 3 


A> A 3 


Maslov(U+,E^) 





-1 


-2 


-3 



17 A non-monotone Maslov index. 



In the LW-SW system and the KdV5 equation, the Maslov index is a monotone function 
of A. (Note however that it is not a monotone function of x.) Here we show an example 
where the Maslov index is not a monotone function of A. It is a slight modification of the 
LW-SW resonance equations. In this case the correlation between the number of roots of 
the Evans function and the value of the Maslov index is no longer apparent. Look at the 
eigenvalue problem 

with L|' h=| '-^-2«M« + 28W» + 2- h 
w xx — cw + 6w(x)w + 2u{x)u 



















L 














w 
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Fig. 11. Longwave-Shortwave problem for the following parameters c = 1, v = 0.2. Top: Maslov 
index as a function of A. Bottom: Evans function as a function of A. The Maslov index is 
A-decreasing. 

with c> 4z/ > 0, 

u(x) = Asech(v / ^ x) and w(x) = 2i/sech 2 (-\/z/a;) 
with A 2 = 2v[c — 4z/), and the requirement c > 4u > 0. 

The spectral problem associated to this operator can be expressed in the form (jl.2p with 
n = 2, 



u(x) 



( u \ 



w 

2u, r 



, and B(x, A) 



A - 2v + 2w{x) 
2u(x) 





2«(:r) 0^ 
-\-c + 6w(x) 








1 



The essential spectrum of L consists of 

] — oo, — c] U [2r/, +oo[. 
The essential spectrum is unbounded from above and below, hence a Morse index cannot 
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Fig. 12. Plot of the Maslov index for the non-monotone example (|17.ip for the case c = 1 and 
v = 0.21. Upper figure shows the Maslov index and the lower figure the Evans function. 

be defined for L. However, we will still be able to compute a Maslov index. The key 
property that leads to non-monotonicity is that the matrix d\B(x, A) is not semi- definite: 
the matrix d\B(x, A) has eigenvalues {0, 0, —1, +1} and so is not semi-definite. 

Results for the case c = 1 and v = 0.21 are tabulated below and shown in Figure [T2l In 
this case there are 5 eigenvalues, but there is no longer a correlation between the Maslov 
index and the number of eigenvalues in a subset of A. 



A 


-c 2u 


D(X) 


+oo + + - + - -oo 


Maslov(A) 


-4 -3 -2 -1 -2 -3 



18 Concluding remarks 

We have only just scratched the surface of the implications of the Maslov index for ho- 
moclinic orbits and solitary waves. Other important questions for Hamiltonian systems 
on four- dimensional phase space are (a) the connection between transversality of the ho- 
moclinic orbit and the parity of the Maslov index, (b) the jump of the Maslov index at 
bifurcations, (c) whether the angles K\ and k 2 in the decomposition n = Ki + k 2 contain 
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other useful information, (d) the role of purely imaginary eigenvalues in the stability of 
solitary waves in KdV5, and (e) the Maslov index of multi-pulse homoclinic orbits. The 
latter question is addressed in the paper [T7j . 

The extension of the Maslov index of homoclinic orbits to phase space with dimension 
greater than four is straightforward in principle but there are some differences in detail. 
First, the dimension of the basic manifolds (A n (K 2n_1 ), MP ^ -1 , G„(M 2n ), A(n), A x (n)) 
jumps a lot when n goes from 2 to 3. 



Manifold 


A 3 (K 6 ) 


MP 19 


G 3 (K 6 ) 


A(3) 


A*(3) 


Dimension 


20 


19 


9 


6 


5 



In this table, RP 19 is the projectification of A 3 (R 6 )- The biggest change in the numerics 
is the difficulty due to the jump in dimension of the Lagrangian Grassmannian. Whereas 
it is 3— dimensional in the case n — 2, it jumps to double that dimension when n — 3. 
The details of the theory and numerics for n > 3 are given in Part 2 [TH] . 
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Appendix 



A Kernel of A^ 2 ) on A 2 (^ 4 ) 

Let A be an arbitrary 4x4 matrix with entries a^. Then, with respect to the standard 
basis (17. 1 \ on A 2 (^ 4 ) the induced matrix is 



A (2) 



an + a-22 


^23 


«24 


-Ol3 


— ai4 





O32 


an + O33 


«34 


012 





— 014 


(Z42 


043 


Oil -|- (Z44 





012 


Ol3 


-031 


021 





022 + 033 


034 


— 024 


— <24i 





«21 


043 


022 + 044 


«23 





— a4i 


3 1 


— Q42 


032 


033 + a 44 



(a-i; 



A constructive proof is given in §2 of [2]. 

Proposition 10 The induced matrix A^ 2 ) satisfies A^uj 
W. if and only if J A is symmetric. 

Proof. An explicit calculation gives 

( a 23 - ai4 ^ 
an + a 33 
043 + 012 

«21 + a 34 
022 + O44 

\-a 41 + a 32 / 



0, where u is defined in 



A (2) 



CO 



On the other hand 



JA 



—a%i — a.32 —033 —034 

— 041 — 042 — 043 — 044 

On 012 031 014 

^21 O22 O23 Q24 
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In order for JA to be symmetric we require 



041 — a 32 j 0,11 — — «33 1 a 12 — — «43 , 0,21 ~ — &34 j «22 — ~ «44 , &23 — °14 • 

These conditions are satisfied if and only if A^uj = 0. ■ 
B The connection between A( 2 ), 3^ and B^ 2 ) 

In this appendix the role of the induced symplectic operator on A 2 (^ 4 ) is explored. 
Using the standard formula for the induced matrix ( 1A-1I) . the induced form of the sym- 
plectic operator is 











-1 


1 





°1 
























1 














-1 




-1 














1 






























1 - 


-1 








Note that J^uj = and uj t J^ = using uj defined in (16. ip . In fact, the kernel of is 
four dimensional Kernel(J (2) ) = span{E 2 , E 5 , E x + E 6 , E 3 + E 4 }, where Ei, . . . , E 6 is the 
standard basis for /\ 2 (IR 4 )- 
Now suppose B is a symmetric matrix and 





A 


= J _1 B - 


= -JB. 






13 + &24 


&34 


644 


-633 


-hi 





-&12 





-614 


b 23 





-&34 


-^22 


-b 2 3 


^13 ~~ ^24 





&23 


&33 


ha 


biA 





&24 — ^13 


-&14 


-644 


bn 





bu 


-^23 





&34 





bi2 


-bn 


^22 


-bi2 


-b 13 - 



The induced matrix A^ 2 ' does not equal the product of the induced matrices for J and B 
but it has the following form 

a ( 2 ) = -j( 2 )b( 2 ) + S, 
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Fig. C.l. The three cases of hyperbolic spectra for constant coefficient Hamiltonian systems on 
R 4 . 

where S is the skew-symmetric matrix 



S = 






bl2 


-&11 


b22 


-6l2 


-6 2 4 - 


-&12 





-&14 


b 2 3 





-b 3 4 


fell 


614 





&24 - &13 


-614 


-644 


-b 2 2 


-^23 


&13 - ^24 





b 23 


^33 


b\2 





&14 


-^23 





&34 


&13 + ^24 


&34 


644 


-^33 


-^34 






The skew-symmetric matrix S has the properties So; = 0, uj t S = 0, and B^ 2 ^ has the 
properties 

B (2) w ^ but B (2) u; is in the kernel of J (2) . 

Hence the property A^u; = is recovered. Moreover, since uj t J^ = 0, we also have 
u T A^ = 0. 

C Hyperbolic subspaces and Lagrangian planes 

For a linear constant-coefficient Hamiltonian system on M 4 , u x = Au, there are three 
cases where the spectrum of A is strictly hyperbolic and they are shown in Figure IC.ll 
The unstable subspace in each case is a Lagrangian plane. This observation is a special 
case of a result of Montaldi [37] . Here we sketch a proof based on eigenvectors. 

Take the first case. The unstable eigenvalues are of the form /i = v ± it with v > and 
r > 0, with eigenvectors £1 ± i^2 where 

Boo(£i + i6) = (j/ + iT)J(£i + i&), 

or 

Boo£i = i/J£i - rJ£ 2 , Boo^ = zaJ6 + rJ£i . 
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Taking the inner product of the first equation with £ 2 and the second with £1 and sub- 
tracting, 

= fjB^fi - gB^z = vgjfr - = 2^ 2 T J6 • 

Now since v > it follows that £jJ£i = and so span{£i, £2} is a Lagrangian subspace. 

A similar calculation verifies the other two cases, with the third requiring the introduction 
of generalized eigenspaces. 

Now, use the fact that Lagrangian subspaces are invariant for the x— dependent system, 
as shown in §U to conclude that hyperbolic sets are Lagrangian manifolds. 



D The induced inner product on A fc (^ 2n ) 



In this appendix the equivalence between the induced inner product [-,-]& on /\ k (M. 2n ), 
which has dimension d, and the standard inner product on M. d is established. Here the 
general case of 2n— dimensional phase space is considered, which will be required in Part 

2 ng. 

With the standard orthonormal basis for M? n , {e l5 . . . , e 2n }, the nonzero and distinct 
members of the set 

(D-l) 



{e^A-'-Ae^ : i u . . . , i k = 1, . . . , 2n } 



form a basis for the vector space A 



k {TU>2n\ 



with exactly d 



distinct elements. 



{n-k)\k\ 

Choose an ordering such as a standard lexical ordering and label the nonzero distinct ele- 
ments in the set (ID- II) by Ei, . . . , E^. Then, any element U G /\ k (R 2n ) can be represented 
as U = J2j=iUj~Ej. The inner product (-, •) on M 2n induces an inner product on each 
vector space /\ k (R 2n ) as follows. Let 



U = ui A ■ ■ ■ A u fe and V = v x A ■ ■ ■ A v fe , Uj, v,- e M. 2n , V i, j — 1 



, k , 



be any decomposable fc-forms. A k— form is decomposable if it can be written as a pure 
form: a wedge product between k linearly independent vectors in M 2 ™. The inner product 
of U and V is defined by 



[U,V] fc :=det 



(Ui, Vi) ■ ■ • (Ui, v fc ) 



(U fc , Vi) ■ • • (life, Vfc) 



U, V e A fc (R 2n ) 



(D-2) 



Since every element in /\ k (R 2n ) is a sum of decomposable elements, this definition extends 
by (bi)-linearity to any A;-form. Using the orthonormality of the induced basis 



|Ei,E 



j\k 



1 if % = j 
OiH^j 



48 



the inner product between two elements U = Ef=i and V = EjLi ^j'Ej is 

[U, V] fc = [[Eti UiE h Eti fe = Eti EU u * v i I E - E ;l* 



EiiEW:=(U,V) 



d : 



yielding the equivalent representation 

[U, V] fe = (U, V), , U, V G A fc (K 2n ) • (D-3) 



E Plus and minus subspaces in A 2 (^ 4 ) 

Consider R 4 with the standard symplectic and volume forms. Let V = span-f^,^} be a 
two-dimensional oriented subspace of M 4 . 

An oriented subspace is defined as follows. A subspace is an equivalence class of bases; 
that is, span{£i,£2} and span{r/i, 772} represent the same subspace if and only if there is 
an invertible 2x2 matrix m such that [^1^2] = folium.. A subspace is oriented if m 
is restricted to have positive determinant. The oriented subspace span{£i,£ 2 } has one of 
three types 

plus subspace if uj A £1 A £ 2 > 

Lagrangian subspace if uj A £1 A £ 2 = 

minus subspace if uj A £1 A £2 < 0, . 
Proposition 11 The sign o/wAU is an invariant of A3. 3\) . 
Use Proposition [2] to conclude that 

uj A U(x) = constant , 

along solutions of (13.31) . 

There is an interesting connection between Krein signature and the above classification 
of oriented subspaces. Krein signature is a sign which is associated with purely imaginary 
eigenvalues (in the linearization about an equilibrium) or Floquet multipliers (in the 
linearization about a periodic orbit). Consider the case of a simple purely imaginary 
eigenvalue iu, v > 0. Its complex eigenvector ( = £1 + 1^2 satisfies 

BC = uaJC. 
The Krein signature is defined as the sign of 

i<jc,c) = 2<j£ 2 ,a>. 

Now use the identity (J£ 2 , v) = [a;, £ 2 A v] 2 , V v e R 2n to obtain 
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i(JC, C)vol = 2(J£ 2 , 6)vol = [«, 6 A 6] 2 vol = -[u, 6 A 6I2V0I = a; A Ci A e 2 ■ 

This observation also emphasizes the fact that a choice of orientation underlies the defi- 
nition of Krein signature. 



F A spectral problem with sech potential 



This appendix establishes the basic properties of the ODE eigenvalue problem 



+ 12 sech xd> = Hi 



(F-l) 



in the set J?T :={kgIR|k>0}. The solutions of this ODE can be determined explicitly. 
The eigenvalues are re = 1, 4, 9. For all k G Jff \ {1,4, 9}, the two functions 



(^{x, re) = e ±v ^ x (±a + ai tanh(x) ± a 2 tanh 2 (x) + tanh 3 
are linearly independent, where 



x 



(F-2) 



\J K 1 

a = — (4 -re), a x = -(2re-3), a 2 = -<\/re. 

The eigenvalues can be verified by explicit calculation. That k = 1,4,9 are the only 
eigenvalues in Jtf, and that are linearly independent is verified by computing the 
Wronskian 

W(x, k) = det 

It is easily verified that W x = and so W (x, re) is independent of x. Evaluate at x = 



W(0, re) = det 



cii + a ^ s /K ai + ao\/re 



2a (ai + a v^) 



Substituting for a and a 1; 



W(0,re) = — ^(re-l)(re-4)(re-9). 



Hence <f> are linearly independent for all k e \ {1, 4, 9}. 
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The eigenfunctions are 

4>(x, 1) = sech(x)(4 — 5sech 2 (x)) , when k — 1 
0(x,4) = tanh(x) sech 2 (x) , when n = 4 
(fi(x, 9) = sech 3 (a;) , when n = 9 . 
modulo an arbitrary multiplicative constant. 

G Attractivity of the Lagrangian Grassmannian A(2) 

One of the advantages of subtracting off the growth rate at infinity in the equations 
on A 2 (K 4 ), as in (j!2.ip . is that the Lagrangian Grassmannian becomes an attracting 
invariant manifold. When A(2) is attractive, one has greater freedom in choosing the 
numerical integration scheme. 

To prove attractivity, consider the integration of the 2— form representing the unstable 
subspace U + (x, A) 

^-U+ = A< 2 >(x,A)U+ U + 6A 2 (1 4 ) ~L<x<+L. 
ax 

Introduce the transformation 

U + 0r,A) =e ff +W x U + (x,A) 

where cx + (A) is the sum of the eigenvalues of Aoo(A) with positive real part. Then U + 
satisfies 

^-U+ = [A (2) (x,A) -a + (A)I]U+ — L < x < +L (G-l) 
The Lagrangian Grassmannian is the set 

U+ A U + = and u> A U + = . 
When evaluated on the differential equation (1G-1I) these invariants satisfy 

f U+ A U+ = f U+ A U+ + U+ A f U+ 

ax ax ax 

= A< 2 )U + A U+ + U+ A A( 2 )U+ - 2a + U+ A U+ 
= Trace (A) U+ A U+ - 2a + U+ A U+ 
= -2a , U+ A U+ , 
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since Trace (A) = 0. A similar calculation with oj A U + yields 
d 



dx 



ijj A U + = oj A (A (2) - a+I)U+ = -A (2) u; A U+ - a+ w A U + = -a+ u A U + , 



using the fact that A^uj = 0, which is proved in Appendix |A], This proves that 



U+AU+fx) = e- 2CT+x U + AU H 



and wAU" 



\x 



l 'cl> AXJ 



for x > — L 



The eigenvalue cx + is real and positive. Hence when integrating the unstable subspace U + 
along the Lagrangian Grassmannian, both U + A U + and oj A U + (x) are exponentially 
attracted to the zero set. Therefore a special integrator is not required for maintaining 
the constraints; a standard Runge-Kutta algorithm is quite satisfactory. 



H The existence of at least one negative eigenvalue 

Consider the linear operator 

Jzf0 := (f) xxxx - P(j) xx + a(x)(j) , (H-l) 

introduced in (114. ip with a(x) = c — (q + l)(f>(x) q and satisfying (114.51) . Assume 

P + 2c > and < c < 1 or P > and c> . (H-2) 

The essential spectrum for this problem is non-negative. Here it is proved that Jzf has at 
least one negative eigenvalue in the point spectrum. 

Multiply (I14.5P by the basic state (f>(x), 



= C(J) 2 



= c<f) 2 - (P + 2c) 4><j) xx + 2c$$ xx + 4>4> xxxx 

= c(4> + 4> xx ) 2 - (P + 2c) 4>4> xx - c4> 2 xx + 4>4> xxxx . 

Hence integrating, using the fact that <fi and its derivatives decay exponentially as x 
±oo, and the hypotheses ( 1H-2D 



<f) g+2 dx = 
or if P > and c> 0, 



l xx ) 2 dx + (P + 2c) 



:dx + (1 - c) 



b 2 

XX 



dx > . (H-3) 



& +2 d 



x 



[c4> 2 + P4> 2 X + gj dx > o . 



(H-4) 
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To prove that (1H-1I) has a negative eigenvalue, we will show that the quadratic form 
(u,Jifu) is negative when u = <p where (u,Jfu) := Jzfudx. Now 



(u, ££u) 



U=<f) 



J— oo ( $ ( ^Pxxxx - P<Pxx + a 0)0)) dx 

Ho(& + ^8 + c ^ 2 ) dx - (g + 1) 0« +2 d* 

-9 Ho 9+2 dx , 



using (E3j) in the last line. It follows from (IhPl or (jHgj) that (0, Jgf0) < 



I Check of hypothesis [H] for the KdV5 system 

In this appendix, the details are given of the proof that A(x, A) = J -1 B(:r, A) for KdV5, 
with B(x, A) defined in (114.21) . satisfies Hypothesis El 

First set 



(l-A)i 

When A is large and negative, s is a small parameter. This parameter will be used to 
obtain series expansions of the eigenvalues and eigenvectors. 

The characteristic polynomial of A oc (A) is 

= det[XI - A oc (A)] = X 4 - PX 2 + — . 

This polynomial is a biquadratic and for s small it has four complex roots, one in each 
quadrant. Let 9(s) be the eigenvalue in the right-upper quadrant. Its Taylor expansion is: 

OU) = _L (l - I ip s 2 _ lp2 4 _ J_ ip 3 s 6 + q f 8 

V ; V 4 32 128 V 



The other eigenvalues are 9(s), —9(s),—9(s). The eigenvector associated with 9(s) is : 

s 4 e 



v is] 




1 

V » 2 I 



A Taylor expansion of this eigenvector is: 

v(s) = svi + isv 2 + O (s 6 ) . 
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with 



2 



and 



v 2 (s 



V 



-2_ l s 2p2 



/ 



(3fo(s), £sv(s)) is a basis of the unstable space. Let V uns t the matrix whose columns are 
(s) and ^sv(s)). 

The eigenvector associated to — 9(s) is: 



w(s) = s 



1 

V 2 / 



(tw(s),3w(s)) is a basis of the unstable space. Let U st the matrix whose columns are 
5Rw(s) and Q : w(s)). 



The matrix ( V unst \U st ) is not a symplectic matrix but V(s) 
—UstiV^^JUst)" 1 , is. Besides, we have: 



V unst \V st , With 



J_ 

2s 



.1 c2) 

4 ' 



32 



2 h -1- 16 a r 



0(s e 



We also have 



J( V)J since is a symplectic matrix. 
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B 



Let 

/o 1 (A 



\o o o oy 

We are now able to evaluate: = O (s 2 ) but also 

d o\ 

1 




Therefore, {V 1 BV)® e x = O (s 2 ). Hence, as R(x, A) = A(x, A) - A^A) = (1 - a(x))B 
and as |1 — a(x)\ < Cie~° 2 ^, this proves that Hypothesis [HI of Proposition is satisfied. 



J Hamiltonian evolution equation for LW-SW equations 



The LW-SW equations (116.21) can be expressed in Hamiltonian form as follows. Let 



K 



\ 
"5° 

o o -a. 



H(Z) = J (u 2 x + vl+ -w 2 x - w(w 2 + u 2 + v 2 ) + -cw 2 + u(u 2 + t; 2 ) ) 

— oo 

with Z = (u, v, w). Then the system becomes 

Z t = KVH(Z) , 
since, with respect to an L 2 (M) inner product, 



VH(Z) 



\H W j 



—2u xx — 2uw + 2vu 
—2v xx — 2vw + 2vv 
—w xx + cw — 3w 2 — u 2 — v 2 
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K The reduced eigenvalue problem associated with LW-SW equations 



The two-dimensional ODE (116.71) that arises in the reduced problem for LW-SW resonance 
can be written in the form 

1 



v xx + 2z/sech {\fvx) v — [u — ~\ \ v . 

ODEs of this type can be solved explicitly as noted in Appendix[F] The essential spectrum 
is the semi-infinite interval a ess (L) = [2u, +00). Now suppose that A < 2v. Then the 
system at infinity is hyperbolic and one can explicitly construct the solutions (v + ,v~) 
which give the solutions for the stable and unstable subspace 



v ± (x 



A) = e^^iTV + tanh(v^x)) , H=\\l 



The Evans function can be obtained from 



D{\) = det 



The Maslov index is 



v + (x; A) v (x; A) 
Maslov(A) = 



2^ 



A 



v 



2v 



x=0 



1 if A < 

if < A < 2v 
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